Magnetism and Thermodynamics of Spin- 1/2 Heisenberg Diamond Chains in a 

Magnetic Field 



O 

o 

(N 

u . 

Ph. 
< 



o 
o 



> 
o 
in 



o 
o 



X 



Bo Gu and Gang Su* 

College of Physical Sciences, Graduate University of Chinese Academy of Sciences, P. O. Box 4588, Beijing 100049, China 

The magnetic and thermodynamic properties of spin- 1/2 Heisenberg diamond chains are investi- 
gated in three different cases: (a) Ji, J2, J3 > (frustrated); (b) Ji, J3 < 0, J2 > (frustrated); 
and (c) Ji, J2 > 0, J3 < (non-frustrated), where the bond couphng Ji {i = 1,2,3) > stands 
for an antiferromagnetic (AF) interaction, and < for a ferromagnetic (F) interaction. The density 
matrix renormalization group (DMRG) technique is invoked to study the properties of the system 
in the ground state, while the transfer matrix renormalization group (TMRG) technique is applied 
to explore the thermodynamic properties. The local magnetic moments, spin correlation functions, 
and static structure factors are discussed in the ground state for the three cases. It is shown that 
the static structure factor S{q) shows peaks at wavevectors q = avr/S (a = 0, 1, 2, 3, 4, 5) for different 
couphngs in a zero magnetic field, which, however in the magnetic fields where the magnetization 
plateau with m = 1/6 pertains, exhibits the peaks only at g = 0, 27r/3 and 47r/3, which are found 
to be coupUngs-independent. The DMRG results of the zero- field static structure factor can be 
nicely fitted by a linear superposition of six modes, where two fitting equations are proposed. It is 
observed that the six modes are closely related to the low-lying excitations of the system. At finite 
temperatures, the magnetization, susceptibihty and specific heat show various behaviors for different 
couphngs. The double-peak structures of the susceptibility and specific heat against temperature 
are obtained, where the peak positions and heights are found to depend on the competition of the 
couplings. It is also uncovered that the XXZ anisotropy of F and AF couplings leads the system of 
case (c) to display quite different behaviors. In addition, the experimental data of the susceptibility, 
specific heat and magnetization for the compound Cu3(C03)2(OH)2 are fairly compared with our 
TMRG resuhs. 

PACS numbers: 75.10.Jm, 75.40.Cx 



I. INTRODUCTION 

Low-dimensional quantum spin systems with compet- 
ing interactions have become an intriguing subject in the 
last decades. Among many achievements in this area, 
the phenomenon of the topological quantization of mag- 
netization has attracted much attention both theoreti- 
cally and experimentally. A general necessary condi- 
tion for the appearance of the magnetization plateaus 
has been proposed by Oshikawa, Yamanaka and Affleck 
(OYA) [l|, stating that for the Heisenberg antiferromag- 
netic (AF) spin chain with a single-ion anisotropy, the 
magnetization curve may have plateaus at which the 
magnetization per site m is topologically quantized by 
n{S — m) = integer, where S is the spin, and n is the 
period of the ground state determined by the explicit 
spatial structure of the Hamiltonian. As one of fasci- 
nating models which potentially possesses the magne- 
tization plateaus, the Heisenberg diamond chain, con- 
sisting of diamond-shaped topological unit along the 
chain, as shown in Fig. [H has also gained much atten- 
tion both experimentally and theoretically (e.g. Refs. 

It has been observed that the compounds, 
A3Cu3(P04)4 with A = Ca, Sr[2, and Bi4Cu3V20i4[3] 
can be nicely modeled by the Heisenberg diamond chain. 
Another spin-1/2 compound Cu3Cl6(H20)2-2H8C4S02 
was initially regarded as a model substance for the 
spin-1/2 diamond chain Q, but a later experimental 



research reveals that this compound should be described 
by a double chain model with very weak bond alterna- 
tions, and the lattice of the compound is found to be 
Cu2Cl4-H8C4S02 0. Recently, Kikuchi et a/.[| have 
reported the experimental results on a spin-1/2 com- 
pound Cu3 (003)2 (0H)2, where local Cu^+ ions with 
spin S* = 1/2 are arranged along the chain direction, and 
the diamond-shaped units consist of a one-dimensional 
(ID) lattice. The 1/3 magnetization plateau and the 
double peaks in the magnetic susceptibility as well as 
the specific heat as functions of temperature have been 
observed experimentally @, 0|, which has been discussed 
in terms of the spin-1/2 Heisenberg diamond chain with 
AF couphngs Ji, J2 and J3 > 0. 

On the theoretical aspect, the frustrated diamond spin 
chain with AF interactions Ji , J2 and J3 > was studied 
by a few groups. The first diamond spin chain was ex- 
plored under a symmetrical condition Ji = J3 Q . Owing 
to the competition of AF interactions, the phase diagram 
in the ground state of the spin-1/2 frustrated diamond 
chain was found to contain different phases, in which the 
magnetization plateaus at m = 1/6 as well as 1/3 are 
predicted d, [13, Another frustrated diamond 
chain with ferromagnetic (F) interactions Ji , J3 < and 
AF interaction J2 > was also investigated theoretically, 
which can be experimentally realized if all angles of the 
exchange coupling bonds are arranged to be around 90°, 
a region where it is usually hard to determine safely the 
coupling constants and even about their signs [3, ITsj . 
Despite of these works, the investigations on the Heisen- 



2 



berg diamond spin chain with various competing inter- 
actions are stih sparse. 



Motivated by the recent experimental observation on 
the azurite compound Cu3(C03)2(OH)2[l, 0], we shall 
explore systematically the magnetic and thermodynamic 
properties of the spin-1/2 Heisenberg diamond chain with 
various competing interactions in a magnetic field, and 
attempt to fit into the experimental observation on the 
azurite in a consistent manner. The density matrix renor- 
malization group (DMRG) as well as the transfer ma- 
trix renormalization group (TMRG) techniques will be 
invoked to study the ground-state properties and ther- 
modynamics of the model under interest, respectively. 
The local magnetic moments, spin correlation functions, 
and static structure factors will be discussed for three 
cases at zero temperature. It is found that the static 
structure factor S{q) shows peaks in zero magnetic field 
at wavevectors q = air/S (a = 0, 1, 2, 3, 4, 5) for different 
couplings, while in the magnetic fields where the mag- 
netization plateau with m — 1/6 remains, the peaks ap- 
pear only at wavevectors g = 0, 27r/3 and 47r/3, which 
are found to be couplings- independent. These informa- 
tion could be useful for further neutron studies. The 
double-peak structures of the susceptibility and specific 
heat against temperature are obtained, where the peak 
positions and heights are found to depend on the com- 
petition of the couplings. It is uncovered that the XXZ 
anisotropy of F and AF couplings leads the system with- 
out frustration (see below) to display quite different be- 
haviors. In addition, the experimental data of the sus- 
ceptibility, specific heat and magnetization for the com- 
pound Cu3 (003)2 (011)2 are fairly compared with our 
TMRG results. 



The rest of this paper is outlined as follows. In Sec. II, 
we shall introduce the model Hamiltonian for the spin- 
1/2 Heisenberg diamond chain with three couplings Ji, 
J2 and J3, where three particular cases are identified. In 
Sec. Ill, the magnetic and thermodynamic properties of 
a frustrated diamond chain with AF interactions Ji, J2 
and J3 > will be discussed. In Sec. IV, the physi- 
cal properties of another frustrated diamond chain with 
F interactions Ji, J3 < and AF interaction J2 > 
will be considered. In Sec. V, the magnetism and ther- 
modynamics of a non-frustrated diamond chain with AF 
interactions Ji, J2 > and F interaction J3 < will 
be explored, and a comparison to the experimental data 
on the azurite compound will be made. Finally, a brief 
summary and discussion will be presented in Sec. VI. 
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FIG. 1: (Color online) Sketch of the Heisenberg diamond 
chain. The bond interactions are denoted by Ji, J2, and J3. 
Three cases will be considered: (a) Ji, J2, J3 > (a frustrated 
diamond chain); (b) Ji, J3 < 0, J2 > (a frustrated diamond 
chain with competing interactions); and (c) Ji, J2 > 0, J3 < 
(a diamond chain without frustration). Note that Ji > 
stands for an antiferromagnetic interaction while Ji < for a 
ferromagnetic interaction, where i = 1,2, 3. 



II. MODEL 

The Hamiltonian of the spin-1/2 Heisenberg diamond 
chain reads 

L/3 

H — ■ 83^-1 + J2S3i_i • S3i 

i=l 

+ Jz^Zi-2 ■ S3i + J3S3i_i • 83^+1 
L 

+JiS3,-S3,+i)-H.^Sj, (1) 

where Sj is the spin operator at the jth site, L is the to- 
tal number of spins in the diamond chain, (j = 1, 2, 3) 
stands for exchange interactions, H is the external mag- 
netic field, g^B = 1 and ks = 1- Ji > represents the 
AF coupling while Jj < the F interaction. There are 
three different cases particularly interesting, as displayed 
in Fig. [U which will be considered in the present paper: 
(a) a frustrated diamond chain with Ji, J2, J3 > 0; (b) 
a frustrated diamond chain with competing interactions 
Jiy J3 < 0, J2 > 0; and (c) a diamond chain without 
frustration with Ji, J2 > 0, J3 < 0. It should be re- 
marked that in the case (c) of this model, since the two 
end points of the J2 bond represent the two different lat- 
tice sites, it is possible that Ji and J3 can be different, 
even their signs. 

The magnetic properties and thermodynamics for 
the aforementioned three spin-1/2 Heisenberg diamond 
chains in the ground states and at finite temperatures 
will be investigated by means of the DMRG and TMRG 
methods, respectively. As the DMRG and TMRG tech- 
niques were detailed in two nice reviews (l6l. [l7| . we 
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shall not repeat the technical details for concise. In the 
ground-state calculations, the total number of spins in 
the diamond chain is taken at least as L = 120. At finite 
temperatures, the thermodynamic properties presented 
below are calculated down to temperature T = 0.05 (in 
units of |Ji|) in the thermodynamic limit. In our calcu- 
lations, the number of kept optimal states is taken as 81; 
the width of the imaginary time slice is taken as e = 0.1; 
the Trotter-Suzuki error is less than lO^'^; and the trun- 
cation error is smaller than 10~^. 



III. A FRUSTRATED HEISENBERG DIAMOND 
CHAIN (Ji, Ja, Js > 0) 

A. Local Magnetic Moment and Spin Correlation 
Function 

Figure [Ha) manifests the magnetization process of a 
frustrated spin- 1/2 Heisenberg diamond chain with the 
couplings satisfying Ji : J2 : J3 = 1 : 2 : 2 at zero tem- 
perature. The plateau of magnetization per site m — 1/6 
is observed. According to the OYA necessary condition 
[l|, the m = 1/6 plateau of a spin-1/2 Heisenberg chain 
corresponds to the period of the ground state n = 3. 
Beyond the magnetization plateau region, the magnetic 
curve goes up quickly with increasing the magnetic field 
H. Above the upper critical field, the magnetic curve 
shows a s-like shape. To further look at how this magneti- 
zation plateau appears, the spatial dependence of the av- 
eraged local magnetic moment (Sj) in the ground states 
under different external fields is presented, as shown in 
Fig. El^b). It is seen that in the absence of external 
field, the expectation value {Sp changes its sign at every 
three sites within a very small range of (—10"'^, 10"'^) 
because of quantum fluctuations, resulting in the magne- 
tization per site m — J2j=i{^j)/^ = 0- i^j) increases 
with increasing the magnetic field, and oscillates with 
increasing j, whose unit of three spins is gradually di- 
vided into a pair and a single, as displayed in Fig. 
At the field H/Ji = 1.5, as demonstrated in Fig. EJd), 
the behavior of {Sp falls into a perfect sequence such as 
{...,{Sa,Sa,Sb),...} with Sa = 0.345 and Sb = -0.190, 
giving rise to the magnetization per site m = 1/6. In ad- 
dition, such a sequence remains with increasing the mag- 
netic field till H/Ji = 2.5, implying that the to = 1/6 
plateau appears in the range of H/Ji = 1.5 ~ 2.5, as 
manifested in Fig. ^a). When the field is promoted fur- 
ther, the sequence changes into a waved succession with 
smaller swing of {Sa — Sb), as shown in Fig. [2l^e), which 
corresponds to the fact that the plateau state of to = 1/6 
is destroyed, and gives rise to a s-like shape of M{H). It 
is noting that when the plateau state of m = 1/6 is de- 
stroyed, the increase of m at first is mainly attributed to 
a rapid lift of Sb, and later, the double Sa start to flimsily 
increase till Sa — Sb — 0.5 at the saturated field. 

The physical picture for the above results could be un- 
derstood as follows. For the to = 1/6 plateau state at 
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FIG. 2: (Color online) For a spin-1/2 frustrated Heisenberg 
diamond chain with fixed couplings Ji : J2 : J3 = 1 : 2 : 2, 
(a) the magnetization per site m as a function of magnetic 
field H in the ground states; and the spatial dependence of 
the averaged local magnetic moment in the ground states 
with external field (b) H/ J\ — 0, (c) 1, (d) 1.5 and 2.5, and 
(e) 4. 



Ji : J2 : J3 = 1 : 2 : 2, we note that if an approximate 
wave function defined bvp^ 



1 



— ^=(2| T3i-2T3i-li3i) ± I T3i-2i3i-lT3i) 

V6 

± I i3»-2T3.-lT3^)), (^-l,...,L/3) 



(2) 



where tj (ij) denotes spin up (down) on site j, is applied, 
one may obtain {tlJi\S§,_2\tpi) = 1/3, {'^P^\S|.^_J^\'^|J^) = 1/3, 
(tpilS^^lipi) = —1/6, giving rise to a sequence {..., (i, |, 
— |), ...}, and TO = (|-f | — ^)/3 = 1/6, which is in agree- 
ment with our DMRG resuhs {..., (0.345, 0.345, -0.190), 
...}. This observation shows that the ground state of this 
plateau state might be described by trimerized states. 

Let _ffci and Hc2 be the lower and upper critical mag- 
netic field at which the magnetization plateau appears 
and is destructed, respectively. For < H < Hc2, the 
magnetization to — nip — 1/6, namely, the system falls 
into the magnetization plateau state. For < H < H^^ 
and Ji : J2 '■ Jz = 1 '■ 2 : 2, the magnetization curve 
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shows the fohowing behavior 
m{H) = ,n,{^)[l+a,{l-^)-a^{l-^f'\ (3) 

J^ci J^ci J^ci 

where Hc^/ Ji — 1.44, the parameters ai = 2/3 and a2 = 
1. Obviously, when H — 0, m = 0; H — H^^^ m = nip. 
A fair comparison of Eq. ([3]) to the DMRG results is 
presented in Fig. [3)^ a). For Hc2 < H < Hg and Ji : J2 '■ 
J3 = 1 : 2 : 2, where Hs is the saturated magnetic field, 
the magnetization curve has the form of 



m{H) = mp + {H-H,,){k, + {Hs-H)[ 
P2 



f3i 



(4) 



where — {nis — mp)/{Ha ~ with the satura- 
tion magnetization, and /3i, /?2 the parameters. One may 



see that when H = Hr. 



H = H, 



rus- 



A nice fitting to the DMRG result gives the parameters 
H.JJi = 3.15, m, = 1/2, H,/Ji = 4.55, /3i = 0.143, 
(32 = 0.178 and = 0.238, as shown in Fig. EJb). 
It should be remarked that from the phenomenological 
Eqs.Q and (|4]), we find that, away from the plateau 
region, the magnetic field dependence of the magnetiza- 
tion of this model differs from those of Haldane-type spin 
chains where m{H) ^ {H — HdY^^- 

To explore further the magnetic properties of the frus- 
trated spin-1/2 Heisenberg diamond chain in the ground 
states with the couplings Ji : J2 ■ J3 = I 2 : 2 at dif- 
ferent external fields, let us look at the static structure 
factor S{q) which is defined as 



Siq) 
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W?" "-"0/ 



(5) 



where q is the wave vector, and (S^Sq) is the spin corre- 
lation function in the ground state. As demonstrated in 
Fig. Hla), in the absence of the external field, S{q) shows 
three peaks: two at g = 7r/3, 57r/3, and one at g = tt, 
which is quite different from that of the spin = 1/2 
Heisenberg AF chain, where S{q) only diverges at g = tt. 
As indicated by Eq. ([5]), the peaks of S{q) reflect the 
periods of the spin correlation function (S^S^), i.e., the 
peaks at 9 = 7r/3 (57r/3) and n reflect the periods of 6 
and 2 for (SjSq), respectively. As shown in Fig. Hfb), 
in the absence of the external fleld, {S^S^) changes sign 
every three sites, which corresponds really to the peri- 
ods of 6 and 2. With increasing the magnetic field, the 
small peak of S{q) at q — tt becomes a round valley while 
the peak at g = 7r/3 (57r/3) continuously shifts towards 
q — 27r/3 (47r/3) with the height enhanced, indicating 
the corruption of the periods of 6 and 2 but the emer- 
gence of the new period G (3,6) for {SjSq), as shown 
in Fig. mjc). At the field H/Ji = 1.5, two peaks shift 
to q = 2tt/3 and 47r/3 respectively, and merge into the 
peaks already existing there, showing the existence of 



.0.10 
0.05 



— Eq.(3) 

m = 1/6, H =1.44, 

p ' Cl 

a.=2/3, a =1 




0.5 


■(b) 

— Eq.(4) 


— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 


0.4 


. m = 1/6, k = 0.238, 

P c 

H^j=3.15, H^=4.55, 
0.143, 0.178 




0.3 




-□-DMRG 


0.2 






□ -□ 





3.0 



3.5 4.0 

H/J, 

1 



4.5 



FIG. 3: (Color online) For a spin-1/2 frustrated Heisenberg 
diamond chain with fixed couplings Ji : J2 : J3 = 1 : 2 : 2, the 
DMRG results of the magnetization per site m as a function 
of magnetic field H away from the plateau state can be fairly 
fitted by Eqs. (O and ©, (a) for < J/ < He-,, and (b) for 
He-, < H < Hs- 



the period 3 for (S^Sq), as clearly displayed in Fig. (H^d). 
The valley and peaks of S{q) keep intact in the plateau 
state at m = 1/6. When the plateau state is destroyed 
at the field H/Ji — 4, the peaks at q — 27r/3 and 47r/3 
are depressed dramatically while the peaks at q — 7r/3, 
57r/3 and tt appear again with very small heights, reveal- 
ing the absence of the period 3 and the slight presence of 
the period 2 and 6, as shown in Fig. IH^e). At the field 
H/Ji = 4.8, all peaks disappear and become zero, ex- 
cept for the peak at g = 0, which is the saturated state. 
Therefore, the static structure factor S{q) shows differ- 
ent characteristics in different magnetic fields [ll]. On the 
other hand, it is known that S{q) also refiects the low- 
lying excitations of the system. It is thus reasonable to 
expect that the low-lying excitations of the frustrated di- 
amond chain will behave differently in different magnetic 
fields. 

To investigate the zero-field static structure factor S{q) 
in the ground state for the frustrated spin-1/2 diamond 
chains with various AF couplings, the four cases with 
Ji = 1, J3 > 0, and J2 = 0.5, 1, 2, and 4 are shown 
in Figs. [5l^a)-[5l^d), respectively. For J2 — 0.5, as shown 
in Fig. [DJa), S{q) displays a sharp peak at q = tt when 
J3 < 0.5, and three peaks at q = 0, 27r/3 and 47r/3 
when J3 > 0.5. It is shown from the ground state phase 
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FIG. 5: (Color online) The zero-field static structure factor 
S{q) in the ground states for the spin-1/2 frustrated Heisen- 
berg diamond chains with length L — 120, Ji = 1, J3 > 
and J2 taken as (a) 0.5, (b) 1, (c) 2 and (d) 4. 



FIG. 4: (Color online) For a spin-1/2 frustrated Heisenberg 
diamond chain with fixed couplings Ji : J2 : J3 = 1 : 2 : 2, 
(a) the static structure factor S{q) in the ground states under 
different external fields; and the spatial dependence of the 
spin correlation function (SjSq) in the ground states under 
external field (b) H/Ji = 0, (c) 1, (d) 1.5 and 2.5, and (e) 4. 



diagram [l^l that the system is in the spin fluid (SF) 
phase when Ji — 1, J2 — 0.5, J3 < 0.5, and enters into 
the ferrimagnetic (FRI) phase when Ji = 1, J2 = 0.5, 
J3 > 0.5. For J2 = 1, as indicated in Fig. ^h), the 
incommensurate peaks exist, such as the case of J3 = 0.8, 
where the system is in the dimerized (D) phase [io|. For 
J2 = 2, as manifested in Fig. [H^c), S{q) has a sharp 
peak at g = TT and two ignorable peaks at 7r/3 and 57r/3 
when J3 < 1; three sharp peaks at q = 0, 2tt/3 and 47r/3 
at J3 = 1; a round valley (J3 = 1.5) or a small peak 
(J3 = 4) at g = TT and two mediate peaks at 7r/3 and 
57r/3 when J3 > 1. The system with Ji = 1 and J1 — I 
is in the D phase when 1 < J3 < 2.8, and in the SF 
phase when J3 < 1 or J3 > 2.8 10]. For J2 = 4 revealed 
in Fig. Eljd), the situations are similar to that of Fig. 
[5jc), but here only the SF phase exists for the system 
with Ji = 1, J2 = 4 and J3 > [l3|. It turns out that 
even in the same phase, such as the SF phase, the zero- 
field static structure factor S{q) could display different 
characteristics for different AF couplings. In fact, we 
note that the exotic peak of Sijq) has been experimentally 
observed in the diamond-typed compoimd Sr3Cu3(P04)4 



If the spin correlation function for the spin-S chain can 
be expressed as {SjSq) = a{—iye^^'^, where a and (3 are 
two parameters, its static structure factor will take the 
form of 



Siq) 



S{S + 1) ajcosq + e-'^) 
3 cos q + cosh /? ' 



(6) 



which can recover exactly the S{q) of the spin-S AKLT 
chain Siq) = S±lj^^£^^M__ Q with a = (5 + 

1)^/3 and (3 = ln(l + 2/S'). Eq.® has apeak at wavevec- 
tor q = TT. By noting that the zero-field static struc- 
ture factor S{q) for the frustrated diamond chains dis- 
plays peaks at wavevectors q = air/S (a — 0, 1, 2, 3, 4, 5) 
for different AF couplings, the spin correlation func- 
tion (SjS^) could be reasonably divided into six modes 
(Slm+iS^) = Q + aze-(6™+0/3 or (SI^+iS^) = ai{6m + 
l)^^ with i = 6m + I and 1 = 1,2, 6, whose contribu- 
tions to the static structure factor should be considered 
separately [21I. Thus, the static structure factor for the 
present systems could be mimicked by a superposition of 
six modes, which leads to 
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FIG. 6: (Color online) The DMRG results of the zero-field 
static structure factor as a function of wavevector for the spin- 
1/2 frustrated Heisenberg diamond chains are fitted: (a) Ji — 
1, J2 = 0.5, J3 = 0.3 by Eq. (O; (b) Ji = 1, J2 = 0.5, J3 = 4 
by Eq. (O; (c) Ji = 1, J2 = 4, J3 = 0.5, and (d) Ji = 1, 
J2 = 4, Ja = 4 by Eq. (O. 



respectively, depending on which phase the system falls 
into, where ai, ci and /3 are couplings-dependent param- 
eters. 

As presented in Fig. [SI the DMRG results of the static 
structure factor as a function of wavevector are fitted 
by Eqs. ([7]) and ^ for the spin- 1/2 frustrated Heisen- 
berg diamond chains with various AF couplings in zero 
magnetic field. It can be found that the characteristic 
peaks can be well fitted by Eqs. and ([H]), with only 
a slightly quantitative deviation, showing that the main 
features of the static structure factor for the present sys- 
tems can be reproduced by a linear superposition of six 
modes. The fitting results give six different modes in 
general, as shown in Fig. [IJb). 

To further understand the above-mentioned behav- 
iors of the zero-field static structure factor, S{q), we 
have applied the Jordan- Wigner (JW) transformation 
to study the low-lying excitations of the spin-1/2 frus- 
trated Heisenberg diamond chain with various AF cou- 
plings (see Appendix A for derivations). It can be seen 
that the zero-field low-lying fermionic excitation e{k) be- 
haves differently for different AF couplings, as shown in 
Figs. [Zla)-(d). Obviously, these low-lying excitations are 
responsible for the DMRG calculated behaviors of S(q), 
where the positions of minimums of e{k) for different AF 
couplings, as indicated by arrows in Figs. [7Ka)-(c), are 
exactly consistent with the locations of the peaks of zero- 
field static structure factor S{q) shown in Figs. [6Ka)-(c), 
respectively, although there is a somewhat deviation for 
Fig. EKd) and Fig. [Sfd). It also shows that the six 



kin 



FIG. 7: (Color online) The zero-field low-lying fermionic ex- 
citation as a function of wavevector for the spin-1/2 frustrated 
Heisenberg diamond chain with (a) Ji = 1, J2 = 0.5, J3 = 0.3, 
(b) Ji = 1, J2 = 0.5, J3 = 4, (c) Ji = 1, J2 = 4, J3 = 0.5, 
and (d) Ji = 1, J2 = 4, J3 = 4. The arrows indicate the 
locations of minimums of e(fc). 



modes suggested by Eqs. ^ and ^ is closely related to 
the low-lying excitations of the system. 



B. Magnetization, Susceptibility and Specific Heat 

The magnetization process for the spin-1/2 frustrated 
diamond chain with Ji = 1, J3 > 0, and J2 = 0.5 and 
2 is shown in Fig. [SJ^a) and ^h), respectively, where 
temperature is fixed as T/Ji = 0.05. It is found that 
the magnetization exhibits different behaviors for differ- 
ent AF couplings: a plateau at m = 1/6 is observed, in 
agreement with the ground state phase diagram [l^, [Q ; 
for J2 = 0.5, as shown in Fig. [SJa), the larger J3 is, the 
larger the width of the plateau at m = 1/6 becomes; for 
J2 = 2, as presented in Fig. [5]^b), the width of the plateau 
at m = 1/6 becomes larger with increasing J3 < 1, and 
then turns smaller with increasing J3 > 1; for J3 < 1, the 
larger J2, the larger the width of the plateau at m = 1/6; 
for J3 = 2, the larger J2, the smaller the width of the 
plateau at m = 1/6. The saturated field is obviously 
promoted with increasing AF J3 and J2. 

Figures [D^c) and [8jd) give the susceptibility x as a 
function of temperature T for the spin-1/2 frustrated 
diamond chain with Ji = 1, J3 > and J2 = 0.5 
and 2, respectively, while the external field is taken as 
H/Ji = 0.01. For J2 = 0.5, as shown in Fig. [^c), the 
low temperature part of xi^) keeps finite when J3 < 0.5, 
and becomes divergent when J3 > 0.5. As clearly man- 
ifested in the inset of Fig. [5^c), J3 = 0.5 is the critical 
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value, which is consistent with the behaviors of static 
structure factor S{q) in Fig. [5]Ja). For J2 = 2, as shown 
in Fig. El^d), an unobvious double-peak structure at low 
temperature is observed at small and large J3 such as 0.2 
and 3. The temperature dependence of the specific heat 
C with Ji — 1, J3 > and J2 = 0.5 and 2 is shown in 
Fig. ^e) and[5I^f), respectively, while the external field 
is fixed as H/Ji = 0.01. For J2 = 0.5, as given in Fig. 
mje), a double-peak structure of C{T) is observed at low 
temperature for small and large J3 such as 0.3 and 1. 
The case with J2 ~ 2 shown in Fig. [5]Jf) exhibits the 
similar characteristics. Thus, the thermodynamics of the 
system demonstrate different behaviors for different AF 
couplings. As manifested in Fig. EJ the low-lying excita- 
tions behave differently for different AF couplings. The 
double-peak structure of the susceptibility as well as the 
specific heat could be attributed to the excited gaps in 
the low-lying excitation spectrum [2^. 



IV. A FRUSTRATED DIAMOND CHAIN WITH 
COMPETING INTERACTIONS (Ji, J3 < 0, J2 > 0) 



A. Local Magnetic Moment and Spin Correlation 
Function 



Figure ini^a) shows the magnetization process of a frus- 
trated spin- 1/2 Hcisenberg diamond chain in the ground 
states with the couplings Ji : J2 : J3 = — 1 : 4 : —0.5. 
The plateau of magnetization per site m 1/6 is clearly 
obtained. To understand the occurrence of the magne- 
tization plateau, the spatial dependence of the averaged 
local magnetic moment (S^) in the ground states at dif- 
ferent external fields is calculated. It is seen that in the 
absence of the magnetic field, as presented in Fig. (Hb), 
the expectation values of (Sj) change sign every three 
sites within a very small range of (—2 x lO^'', 2 x 10^^), 
resulting in the magnetization per site m = 0. With 
increasing the field, the expectation values of {S^) in- 
crease, whose unit of three spins is gradually divided into 
a pair and a single, as displayed in Fig. [He). At the field 
_ff/| Ji| = 0.05, as illustrated in Fig. [DJd), the behavior of 
shows a perfect sequence such as {..., {Sa, Sb, Sb), ■■■} 
with Sa = 0.496 and Sb — 0.002, giving rise to the mag- 
netization per site m = 1/6. In addition, the sequence 
is fixed with increasing the field until H/\Ji\ =3.2, cor- 
responding to the plateau of m = 1/6. As the field is 
enhanced further, the double Sb begin to rise, and the 
sequence becomes a waved series with smaller swing of 
{Sa — Sb) as revealed in Fig. (H^e), which corresponds to 
the plateau state at m = 1/6 that is destroyed. It is 
noting that the increase of m is mainly attributed to the 
promotion of double Sb, as Sa is already saturated until 
Sa = Sb ~ 0.5 at the saturated field. 

As discussed above, the physical picture of the m = 1/6 
plateau state at Ji : J2 : >/3 = — 1 : 4 : —0.5 can be 
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FIG. 8: (Color online) For the spin-1/2 frustrated Heisenberg 
diamond chains with Ji = 1 and J3 > 0, the magnetization 
process m{H) at temperature T/Ji = 0.05 with (a) J2 = 0.5 
and (b) J2 = 2; the susceptibility xiT) at field H/Ji = 0.01 
with (c) J2 = 0.5 and (d) J2 = 2; the specific heat C{T) at 
field H/Ji = 0.01 with (e) J2 = 0.5 and (f) J2 = 2. 

understood by the following approximate wave function 



T3i-2T3i-li3i) ± I T3i-2i3'i-lT3i))- 



By use of this wave function, we have {4'i\S^i_^2\'^i) — 
1/2, (V'^|5f._l|V'^) = 0, {^^\S|,\iJ^) = 0, leading to a 
sequence of {...,(1/2,0,0),...}, and m = (1/2 + + 
0)/3 = 1/6. This is in agreement with our DMRG re- 
sults {..., (0.496, 0.002, 0.002), ...}. 

The static structure factor S{q) of the frustrated spin- 
1/2 Heisenberg diamond chain with the competing cou- 
plings Ji : J2 : J3 = — 1 : 4 : —0.5 in the ground states 
is considered in different external fields. At zero field, 
shown in Fig. [TUT a). S{q) has three peaks a.t q = 7r/3, 
57r/3 and tt with mediate heights, which reflects the pe- 
riods of 6 and 2 for (SjSq), respectively. As shown in 
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FIG. 9: (Color online) For a spin-1/2 frustrated Heisenberg 
diamond chain with fixed couplings Ji : J2 : J3 = — 1 : 4 : 
—0.5, (a) the magnetization per site m as a function of mag- 
netic field H in the ground states; and the spatial dependence 
of the averaged local magnetic moment (Sj) in the ground 
states with external field (b) H/\Ji\ =0, (c) 0.025, (d) 0.05 
and 3.2, and (e) 3.27. 



FIG. 10: (Color online) For a spin-1/2 frustrated Heisenberg 
diamond chain with fixed couplings Ji : J2 : J3 = — 1 : 4 : 
—0.5, (a) the static structure factor S{q) in the ground states 
with different external fields; and the spatial dependence of 
the spin correlation function {SjSq) in the ground states with 
external field (b) H/\Ji \ = 0, (c) 0.025, (d) 0.05 and 3.2, and 
(e) 3.27. 



Fig. [TUTb). in the absence of the external field, {SjSq) 
changes sign every three sites, corresponding to the pe- 
riods of 6 and 2. With increasing the field, the peak at 
q = TT becomes flat with height depressed forwardly, while 
the peak at q — 7r/3 (57r/3) is divided into two peaks 
shifting oppositely from q = n/S (57r/3) with height de- 
creased, indicating the corruption of the periods of 6 and 
2 and the emergence of new periods for {SjS§), as shown 
in Fig. Mic). At the field H/\Ji\ = 0.05, two shift- 
ing peaks have respectively reached q — 27r/3 and 47r/3, 
and are merged with the existing peaks, showing the oc- 
currence of period 3 for (SjSq), as clearly displayed in 
Fig. fTUTd). The flat and peaks keep constant during the 
plateau state at m = 1/6. When the plateau state is 
destroyed at field H/Ji — 3.27, the peaks at g = 27r/3 
and 47r/3 are depressed sharply, revealing the decay of 
the period 3, as shown in Fig. fTUfe). At the saturated 
field of H/Ji = 3.3, all peaks disappear and become flat 
with the value zero, which is the saturated state. So, the 
static structure factor S{q) shows different characteristics 
in different magnetic fields. Similar to the discussions in 



Fig. m the low-lying excitations of this frustrated di- 
amond chain would also behave differently in different 
magnetic fields. 

To further investigate the zero-field static structure 
factor S{q) in the ground state for the frustrated spin-1/2 
diamond chains with various Ji, J3 < and J2 > 0, two 
cases with Ji = — 1, J3 < and J2 = 1 and 4 are illus- 
trated in Fig. fTTT a) andfTlTb). respectively. For J2 = 1, 
S{q) shows a round peak at g = tt, two sharp peaks at n/3 
and 57r/3 when jJaj < 1, and a very sharp peak at g = 
when I J3I > 1. For J2 = 4, S{q) shows three peaks at 
q = 7r/3, TT and 57r/3 as | J3I < 1, and a very sharp peak 
at q — TT and nearly ignorable peaks at tt/3 and 57r/3 
as jJal > 1. In general, the zero-field static structure 
factor S{q) shows different characteristics with different 
competing couplings, whose exotic characteristics could 
be experimentally observed in the related diamond-typed 
compounds. 

As shown in Fig. [121 the DMRG results of the static 
structure factor as a function of wavevector are fitted by 
Eq.([7]) for the spin-1/2 frustrated Heisenberg diamond 
chains with Ji, J3 < and J2 > 0. It can be found 
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FIG. 11: (Color online) The zero-field static structure factor 
S{q) in the ground states for the spin-1/2 frustrated Heisen- 
gberg diamond chains with length L = 120, Ji = — 1, J3 < 
and J2 taken as (a) 1; (b) 4. 



FIG. 12: (Color onhne) The DMRG results of the zero-field 
static structure factor as a function of wavevector are fitted 
by Eq. ((T} for the spin-1/2 frustrated Heisenberg diamond 
chains with (a) Ji — —1, J2 — 1, J3 = —0.1, and (b) Ji — — 1, 
J2 = 4, J3 = -0.5. 



that the characteristic behaviors can be nicely fitted by 
Eq-O, with only a slightly quantitative deviation, show- 
ing that the main features of the static structure factor 
for the present systems can be captured by a superposi- 
tion of six modes. It is consistent with the fact that the 
spin correlation function for the present systems has six 
different modes, as manifested in Fig. [TUTb). 

Similar to the case (a) with all AF couplings in the 
last section, the characteristics of zero-field S{q) for the 
spin-1/2 frustrated Heisenberg diamond chain with Ji, 
J3 < and J2 > can be further undersood in terms 
of the low-lying excitations of the system (see Appendix 
A). By means of the JW transformation, the zero-field 
low-lying fermionic excitation e{k) of the present case is 
calculated, as shown in Figs. [T5f a)-(b). It is found that 
the zero-field low-lying fermionic excitation e{k) differs 
for different couplings, but the positions of minimums 
of e{k) for different couplings, as indicated by arrows in 
Figs. Unta) and (b), appear to be the same. One may see 
that these positions coincide exactly with the locations of 
peaks of the zero-field static structure factor S{q) mani- 
fested in Figs. [TWa) and (b), respectively, showing that 
our fitting equation is qualitatively consistent with the 
low-lying excitations of the system. 
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FIG. 13: (Color online) The zero-field low-lying fermionic 
excitation as a function of wavevector for the spin-1/2 frus- 
trated Heisenberg diamond chain with (a) Ji = — 1, J2 = 1, 
J3 = -0.1, and (b) Ji = -1, J2 = 4, J3 = -0.5. The arrows 
indicate the locations of minimums of e(fc). 
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B. Magnetization, Susceptibility and Specific Heat 



Figures [MT a) and [T4jb) show the magnetization pro- 
cess for the spin- 1/2 frustrated diamond chain at a finite 
temperature T/\Ji\ = 0.05 with Ji = —1, J3 < 0, and 
J2 = 1 and 4, respectively. It is shown that the magneti- 
zation exhibits different behaviors for different Ji , J3 < 
and J2 > 0. A plateau at m = 1/6 is observed at small 
I J3I; with a fixed J2, the larger | J3I, the smaller the width 
of the plateau at m = 1/6, and the plateau disappears 
when I J3I exceeds the critical value; for a fixed J3, the 
larger J2, the wider the width of the plateau at m = 1/6; 
the saturation field is obviously depressed with the in- 
crease of IJ3I at a fixed J2, and is enhanced with the 
increase of J2 at a fixed J3. 

Figures [TWc) and \\Md) manifest the susceptibility x 
as a function of temperature T for the spin-1/2 frus- 
trated diamond chain with Ji = — 1, J3 < and J2 = 1 
and 4, respectively, where the external field is taken as 
H/\Ji\ = 0.01. For J2 = 1, the low temperature part of 
x(T) keeps finite when IJ3I < 1, and becomes divergent 
when IJ3I > 1. As clearly revealed in the inset of Fig. 
[T4jc), J3 = — 1 is the critical value, which is in agreement 
with the behaviors of static structure factor S{q) shown 
in Fig. [TlTa). For J2 = 4, a clear double-peak structure 
of x{T) is obtained at | J3I =8. The temperature depen- 
dence of the specific heat C with Ji = — 1, J3 < and 
J2 = 1 and 4 is shown in Figs. [Hf e) and [nf f). respec- 
tively, where the external field is fixed as i//| Ji| = 0.01. 
For J2 = 1, a double-peak structure of C(T) is observed 
for the case of | J3I = 0.5. The case with J2 = 4 shown in 
Fig. [HTfl exhibits the similar characteristics. It is also 
found that, owing to the competitions among Ji, J3 and 
J2, the thermodynamics demonstrate rich behaviors at 
different couplings. As reflected in Fig. [Til the low-lying 
excitations behave differently with various F interactions 
Ji, J3 and AF interaction J2, while the excitation gaps 
could induce the double-peak structure in the suscepti- 
bihty as well as in the specific heat [2^. 
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FIG. 14: (Color online) For the spin-1/2 frustrated Heisen- 
berg diamond chains with Ji = — 1 and J3 < 0, the mag- 
netization process m{H) at temperature T/Ji — 0.05 with 
(a) J2 = 1 and (b) J2 = 4; the susceptibility x{T) at field 
H/Ji = 0.01 with (c) J2 = 1 and (d) J2 — 4; the specific heat 
C(r) at field H/Ji = 0.01 with (e) J2 = 1 and (f) J2 = 4. 



V. A DIAMOND CHAIN WITHOUT 
FRUSTRATION (Ji, J2 > 0, J3 < 0) 

A. Local Magnetic Moment and Spin Correlation 
Function 

Figure I15r a) shows the magnetization process of a 
non-frustrated spin-1/2 Heisenberg diamond chain in the 
ground states with the couplings satisfying Ji : J2 '■ J3 = 
1 : 2 : —0.5. The plateau of magnetization per site 
771 = 1/6 is observed. The appearance of the magneti- 
zation plateau can be understood from the spatial de- 
pendence of the averaged local magnetic moment (Sj) in 
the ground states under different external fields. In ab- 
sence of the external field, the expectation values of (S^) 
change sign every one site with a waved swing within a 
very small range of (—2 x 10~*, 2 x 10""*), giving rise to 



the magnetization per site m = 0. Under a finite field, 
every three successive spins have gradually cooperated 
into a pair and a single, as shown in Fig. [TSTc). At 
H/Ji = 0.9, as given in Fig. [TSfd). (Sj) shows a perfect 
sequence such as {...,{Sa, St, St), ■■■} with Sa — 0.393 
and Sb = 0.053, resulting in the magnetization per site 
TO = 1/6. Moreover, the sequence is fixed with the field 
increased until H/ Ji = 1.8, corresponding to the plateau 
state of TO = 1/6. As the field is increased further, double 
Sb begin to increase, and the sequence becomes a waved 
succession with a smaller swing of [Sa — Sb), as mani- 
fested in Fig. fTSTe). which corresponds to the fact that 
the plateau state at m = 1/6 is destroyed. It is observed 
that the increase of m at first is mainly attributed to the 
speedy boost of double Sb, and later, Sa starts to increase 
weakly until Sa = Sb = 0.5 at the saturated field. 
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FIG. 15: (Color online) For a spin-1/2 non-frustrated Heisen- 
berg diamond chain with fixed couplings Ji : J2 : J3 = 1 : 
2 : —0.5, (a) the magnetization per site m as a function of 
magnetic field H in the ground states; and the spatial de- 
pendence of the averaged local magnetic moment (Sj) in the 
ground states with external field (b) H/Ji = 0, (c) 0.7, (d) 
0.9 and 1.8, and (e) 2.5. 



For this non-frustrated case with couplings Ji : J2 '■ 
J3 = 1 : 2 : —0.5, the obtained perfect sequence of 
{...,(0.393,0.053,0.053),...} for the m = 1/6 plateau 
state could be understood by the following approximate 
trimerized wave function 

■01 — 2(21 T3i-2T3j-li3i) ± 2| t3i-2i3i-lT3j) 
± I i3j-2T3i-lT3i))- 



(9) 



According to this function, we find (V'i|'S'3i_2lV'i) = ^/I^j 
mSl,_^^^) = 1/18, (V.I^IJ^.) = 1/18, giving rise to a 
sequence of {..., (7/18, 1/18, 1/18), ...}. It turns out that 
m = (7/18 + 1/18 -f l/18)/3 = 1/6. This observation 
implies that the ground state of the plateau state can 
also be described by the trimerized states. 

The static structure factor S{q) of the non-frustrated 
spin-1/2 Heisenberg diamond chain in the ground states 
with the couplings Ji : J2 : J3 = 1 : 2 : — 0.5 is probed 
in different external fields. As shown in Fig. [rBFa). in 
absence of the external field, S{q) shows a sharp peak at 
q — TT, similar to the behaviors of the S* = 1/2 Heisen- 
berg AF chain, which reflects the period of 2 for {SjSq). 
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FIG. 16: (Color online) For a spin-1/2 non-frustrated Heisen- 
berg diamond chain with fixed couplings Ji : J2 : J3 = 1 : 2 : 
—0.5, (a) the static structure factor S{q) in the ground states 
with different external fields; and the spatial dependence of 
the spin correlation function {SjSq) in the ground states with 
external field (b) H/Ji = 0, (c) 0.7, (d) 0.9 and 1.8, and (e) 
2.5. 



As displayed in Fig. fTHTb). at zero external field, (S'JS'o) 
changes sign every one lattice site, corresponding to the 
period of 2. When the field is increased, the peak at 
q = TT becomes a flat with the height depressed greatly, 
while two new peaks with small heights at q = 2tt/3 and 
47r/3 appear, indicating the corruption of the period of 2 
and the emergence of the new period of 3 for (S^Sq), as 
demonstrated in Fig. \M.c). At the field H/Ji = 0.9, 
the peaks at q = 0, 27r/3 and 47r/3 become sharper. 
The flat and peaks of S{q) keep unchanged during the 
plateau state at m = 1/6. When the plateau state is de- 
stroyed at the field H/Ji = 2.5, the peaks at q = 27r/3 
and 47r/3 are suppressed dramatically, revealing the de- 
cay of the period 3, as shown in Fig. [TBTe). At the field 
H/Ji = 2.8, all peaks disappear, and become a flat with 
the value zero, except for the peak at g = 0, which is the 
saturated state. It can be stated that the static struc- 
ture factor S{q) shows various characteristics in different 
magnetic fields. Similar to what discussed in Figs. [4] 
and [ini the low-lying excitations of this non-frustrated 
diamond chain would also behave differently in different 
magnetic fields. 

The zero- field static structure factor S{q) in the ground 
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state for the present system displays a peak at q = ir with 
different couphngs, but whose static correlation function 
{SjSq) varies with the couphngs. As illustrated in Fig. 
[T71 only are the values of (SjSq) larger than zero pre- 
sented for convenience. In order to gain deep insight 
into physics, for a comparison we also include the static 
correlation function for the S" = 1/2 Heisenberg antifer- 
romagnetic (HAF) chain in Fig^T7]^a), whose asymptotic 
behavior has the form of [H, [24 1 



(S,-So)«(-l)-'' 



(27r)3/2 



J 



(10) 



Eq. ([T0| is depicted as solid lines in Figs. \T7\ To take the 
finite-size effect into account, the length of the diamond 
chain is taken as L = 90, 120 and 160, respectively. As 
revealed in Fig. HTJa), the DMRG result of the 5=1/2 
HAF chain with an infinite length agrees well with the 
solid line. Fig. [TWb) shows that the static correlation 
function for the spin- 1/2 diamond chain with frustrated 
couplings Ji : J2 J3 = 1 : 1-2 : 0.5 decays faster than 
that of the HAF chain. Compared with Figs. [TWc)-(f). 
it can be found that all the static correlation functions 
for the spin- 1/2 non- frustrated diamond chain with AF 
interactions Ji, J2 and F interaction J3 fall more slowly 
than that of the 5=1/2 HAF chain; for fixed AF inter- 
actions Ji and J2, the static correlation functions drop 
more leisurely with increasing the F interaction | J3I; for 
fixed AF interaction Ji and F interaction J3, the static 
correlation functions decrease more rapidly with increas- 
ing the AF interaction J2. 



B. Magnetization, Susceptibility and Specific Heat 

Figures [TST a) and [TSf bl show the magnetization pro- 
cess for the spin- 1/2 non- frustrated diamond chain at a 
finite temperature T/Ji = 0.05 with Ji = 1, J3 < 0, and 
J2 = 0.5 and 2, respectively. It is found that the magneti- 
zation behaves differently with different AF interactions 
Ji, J2 and F interaction J3. A plateau at m = 1/6 is ob- 
tained at small IJ3I; for fixed Ji and J2, the larger IJ3I, 
the narrower the width of the plateau at m = 1/6, and 
after | J3I exceeds a critical value, the plateau aXm— 1/6 
is eventually smeared out; for fixed Ji and J3, the larger 
J2, the wider the width of the plateau at m = 1/6; the 
saturated field is obviously unchanged with changing the 
F interaction J3. The coupling-dependence of the spin- 
1/2 non- frustrated diamond chain with AF interactions 
Ji , J2 and F interaction J3 is similar to that of trimerized 
F-F-AF chains [25]. 

Figures fTSTc) and [TSFd) present the susceptibility x 
as a function of temperature T for the spin-1/2 frus- 
trated diamond chain with Ji = 1, J3 < and J2 = 0.5 
and 2, respectively, where the external field is taken as 
H/Ji — 0.01. A double-peak structure of xC^) is ob- 
served at small F interaction | J3 1 and disappears at large 
I J3 1 . The temperature dependence of the specific heat G 
with Ji = 1, J3 < and J2 = 0.5 and 2 is shown in Figs. 
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FIG. 17: (Color online) The zero-field static correlation func- 
tion (SjSq) versus site j in the ground state for a spin-1/2 
diamond chain with different lengths and various couplings. 
The couplings ratio Ji : J2 : J3 is taken as (a) 1 : 1 : (HAF), 
(b) 1 : 1.2 : 0.5 (frustrated), (c) 1 : 2 : -0.5, (d) 1 : 2 : -2, (e) 
1 : 0.5 : -0.5, (f) 1 : 0.5 : -2. The length is taken as i = 90, 
120 and 160, respectively. 



HUe) andlTHf), respectively, where H/Ji = 0.01. It is 
seen that, when J2 is small, C(T) exhibits only a single 
peak; when J2 is large, a double-peak structure of C{T) 
is observed. In the latter case, the double-peak structure 
is more obvious for small F interaction | J3I, and tend to 
disappear at large IJ3I. Therefore, the thermodynamics 
demonstrate various behaviors with different AF interac- 
tions Ji, J2 and F interaction J3. 



C. Effect of Anisotropy of Bond Interactions 

Some magnetic materials show different behaviors un- 
der longitudinal and transverse magnetic fields, showing 
that the anisotropy plays an important role in the physi- 
cal properties of the system. First, let us investigate the 
XXZ anisotropy of the AF interaction J2 on the proper- 
ties of the spin-1/2 non-frustrated diamond chain with 
the couplings Ji : J2Z : J3 = 1 : 2 : —0.5 for vari- 
ous anisotropy parameter defined by 72 = J2x/J2z — 
J2y/J2z, where the z axis is presumed to be perpendic- 
ular to the chain direction. For 72 > 1, the magneti- 
zation m{H), susceptibility x{T) and specific heat C(T) 
are presented in Figs. [TW a). (b) and (c), respectively. 
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on the thermodynamic properties in comparison to what 
we discussed above. 

Now let us discuss the effect of the XXZ anisotropy of 
J3 < on the magnetic and thermodynamic properties 
of the spin-1/2 non-frustrated diamond chain with the 
couphngs Ji : J2 : Jsz = 1:2: —0.5. Recah that as 
the J2 bond connects two different lattice sites, as shown 
in Fig. [U Ji and J3 can be different, even their signs. 
Define a parameter 73 to characterize the anisotropy as 
73 = Jsx/Jsz = J?,yl Jzz, where the z axis is perpendicu- 
lar to the chain direction. For 73 > 1, the magnetization 
m{H), susceptibility x(T) and specific heat C{T) are de- 
picted in Figs. [TWd). (e) and (f), respectively. With 
increasing 73, it is seen that the width of the plateau at 
TO = 1/6 for along the z direction becomes slightly 
wider, while it goes smaller for H along the x direc- 
tion; the saturation field is not changed with 73 along 
both directions; the peak of x{T) for H along the z di- 
rection at lower temperature side is promoted, and the 
second round peak at higher temperature side is slightly 
depressed, while the situations along the x direction are 
just reverse, namely, the peak at lower temperature side 
is depressed, and the second peak at higher tempera- 
ture side is slightly promoted; the peak of C(T) for H 
along the z direction at lower temperature side leaves al- 
most unchanged, and the second round peak at higher 
temperature side moves slightly to the higher tempera- 
ture side, while C{T) along the x direction coincides with 
that along the z direction. For < 73 < 1, the situation 
just becomes reverse in comparison to what we discussed 
above. 



FIG. 18: (Color online) For the spin-1/2 non-frustrated 
Heisenberg diamond chains with Ji = 1 and J3 < 0, the mag- 
netization process m(H) at temperature T/Ji = 0.05 with 
(a) J2 — 0.5 and (b) J2 — 2; the susceptibility x{T) at field 
H/Ji = 0.01 with (c) J2 = 0.5 and (d) J2 = 2; the specific 
heat C{T) at field H/Ji = 0.01 with (e) J2 = 0.5 and (f) 
J2 = 2. 



With increasing 72, it is found that when the magnetic 
field H is along the z direction, the width of the mag- 
netization plateau at to = 1/6 as well as the saturation 
field are enlarged, while those are more increased for H 
along the x direction than along the z direction; the peak 
of the susceptibility x{T) for H along the z direction 
at lower temperature side is promoted, and the second 
round peak at high temperature side is depressed with a 
little shift, while x{T) for H along the x direction shows 
the similar varying trend; the peak of the specific heat 
C{T) for H along the z direction at lower temperature 
side leaves almost unchanged, and the second round peak 
at high temperature side moves towards the higher tem- 
perature side, while C(T) for H along the x direction 
coincides with those for H along the z direction. For 
< 72 < 1, the anisotropy just shows very reverse effect 



D. Comparison to Experimental Results 

Recently, Kikuchi et al. [6|] have performed a nice 
measurement on a spin-1/2 diamond-chain compound 
0113(003)2(011)2, i.e., azurite. They have observed the 
1/3 magnetization plateau, unambiguously confirming 
the previous theoretical prediction. The two broad peaks 
both in the magnetic susceptibility and the specific heat 
are observed. We note that in Ref. the experimen- 
tal data at finite temperatures are fitted by the zero- 
temperature theoretical results obtained by the exact di- 
agonalization and DMRG methods, while the result of 
the high temperature series expansion fails to fit the 
low-temperature behavior of the susceptibility. In ac- 
cordance with our preceding discussions, by using the 
TMRG method, we have attempted to re-analyse the ex- 
perimental data presented in Ref. to fit the experi- 
ments for the whole available temperature region. 

Our fitting results for the temperature dependence of 
the susceptibility x of the compound Ou3(003)2(OH)2 
are presented in Fig. I20f a). For a comparison, we 
have also included the TMRG result calculated by us- 
ing the parameters given in Ref. Obviously, our 
TMRG results with Ji : J2 : Jsz = I ■ 1-9 : -0.3 and 
JzxjJzz = JzyjJzz = 1-7 fit very well the experimental 
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FIG. 19: (Color online) For the spin-1/2 non-frustrated dia- 
mond chain with the couplings satisfying Ji : J2z : J3 = 1 : 
2 : —0.5 for various anisotropy 72 = JixjJ^z > 1: (a) the 
magnetization process m{H) at temperature T/Ji = 0.05; 
(b) the susceptibility x{T at field H/Ji = 0.01; and (c) the 
specific heat C{T) at field H/Ji = 0.01. For the spin-1/2 
non-frustrated diamond chain with couplings Ji : J2 : Jsz ~ 
1:2: —0.5 for various anisotropy 73 — Jsx/Jsz > 1: (d) the 
magnetization process m{H) at temperature T/Ji — 0.05; 
(e) the susceptibility x{T) at field H/,h = 0.01; and (f) the 
specific heat C(r) at field H/Ji^ 0.01. 



data of Xi s-^d the two round peaks at low temperatures 
are nicely reproduced, while the result with Ji : J2 : 
J3 = 1 : 1.25 : 0.45 obtained in Ref. cannot fit the 
low-temperature behavior of x (26j . On the other hand, 
the fitting results for the temperature dependence of the 
specific heat C{T) of the compound Cu3(C03)2(OH)2 
are shown in Fig. I^Uf b). The lattice contribution, which 
is included in the raw experimental data in Ref. [6|, is 
subtracted according to C(T) = CsxpiT) — aT^, where 
a is a parameter. Obviously, our TMRG result with the 
same set of parameters Ji : J2 : J32 = 1 : 1.9 : —0.3 
and J^x/Jzz = Jzyl Jzz = 1-7 fits also remarkably well 



the experimental data of C{T), and the two round peaks 
at low temperatures are nicely reproduced, while the re- 
sult with Ji : J2 : J3 = 1 : 1.25 : 0.45 given in Ref. ^ 
cannot fit the low-temperature behavior of C(T), even 
qualitatively. In addition, the sharp peak of C(T) ex- 
perimentally observed at temperature around 2K cannot 
be reproduced by both sets of the coupling parameters, 
which might be a three-dimensional long-range ordering 
due to interchain interactions. The fitting results for the 
magnetization m{H) of the compound Cu3(C03)2(OH)2 
are shown in Fig. I^OTc). We would like to point out 
that the quantitative fitting by our above parameters to 
the width of the plateau is not so good, but the quali- 
tative behavior is quite consistent with the experiments 
both in the transverse and longitudinal magnetic fields, 

say, H^}^ > H^, ijjj < ^ci^ ^^'^ the saturation field is 
fixed along both directions, suggesting that our fitting 
parameters capture the main characteristics. It is worth 
pointing out that if the anisotropy ratio is increased up 
to Jsx/Jsz — Jsy/Jsz — 2.5 with the same couplings 
Ji ■■ J2 ■■ Jsz = I ■■ 1-9 : -0.3, the width of the 1/3 
plateau ior H \\ b will be decreased to about one-half of 
that for H ±b. 

Therefore, our calculations show that (i) the best cou- 
plings obtained by fitting the experimental data of the 
susceptibility for the azurite could be Ji : J2 : Jzz — 1 : 
1.9 : —0.3 with the anisotropic ratio for the ferromag- 
netic interaction J^x/Jzz = JzylJzz = 1-7, where z _L 6; 

(ii) the compound may not be a spin frustrated magnet; 

(iii) the double peaks of the susceptibility and the specific 
heat are not caused by the spin frustration effect, but by 
the two kind of gapless and gapful excitations owing to 
the competition of the AF and F interactions. 

One might argue that for this diamond chain com- 
pound, from the point of the lattice distance it is unlikely 
that J\ is AF without XXZ anisotropy while J3 is F with 
strong XXZ anisotropy. We may offer another possibil- 
ity to support our findings, namely, the case of J\ and 
J3 with opposite signs is not excluded from the lattice 
structure of the compound. A linear relationship exists 
between the exchange energy and the metal-ligand-metal 
bridge angle: the coupling energy, positive (ferromag- 
netic) at angles near 90°, becomes increasingly smaller 
(more antiferromagnetic) as the angle increases [23|. As 
the ferromagnetic coupling J3 is determined by fitting 
the experimental low-temperature behaviors of x(T) and 
C (T) , this fitting coupling parameters should not be im- 
possible if one considers the angle of Ji bridge to keep the 
antiferromagnetic coupling while the angle of J3 bridge 
to induce the ferromagnetic coupling. On the other 
hand, we note that there is another compound with Cu 
ions, Cu2(abpt)(S04)2(H20)-H20, whose g factors in XY 
plane are different from that in z direction [28]. Besides, 
someone might argue that the condition J2 ^ Ji , | J3 1 is 
necessary to explain the double peak behavior of the dia- 
mond chain. In fact, such an argument is not necessarily 
true, as manifested in Fig. [T5f c). where the double peaks 
of x(T) at low temperatures can also be produced with 
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FIG. 20: (Color online) A comparison of experimental re- 
sults for (a) the magnetic susceptibility, (b) the specific heat 
and (c) the magnetization process for the spin- 1/2 diamond 
compound Cu3(C03)2(OH)2 with the TMRG results. The 
experimental data are taken from Ref. @). See the context 
for details. 



the parameters Ji : J2 ■ J3 — 1 ■ 0.5 : —0.1. In other 
words, the double-peak behavior of the diamond chain 
may not depend on whether J2 ^ Ji, J3 or not, but may 
be strongly dependent on the competition of AF and F 
interactions, as discussed above. 



VI. SUMMARY AND DISCUSSION 

In this paper, we have numerically studied the mag- 
netic and thermodynamic properties of spin- 1/2 Heisen- 
berg diamond chains with three different cases (a) Ji, 
-^2, -^^3 > (frustrated), (b) Ji, J3 < 0, J2 > (frus- 
trated), and (c) Ji, J2 > 0, J3 < (non- frustrated) by 
means of the DMRG and TMRG methods. In the ground 
states, the local magnetic moment, spin correlation func- 
tion, and static structure factor are explored. The static 
structure factor S{q) at zero field shows peaks at wave 



vector q = 0, 7r/3, 2tt/3, tt, 47r/3 and 57r/3 for different 
couplings, in which the peaks at g = 0, 27r/3 and 47r/3 
in the magnetization plateau state with to = 1/6 are ob- 
served to be couplings independent. The DMRG results 
of the zero-field static structure factor can be nicely fitted 
by a linear superposition of six modes, where two fitting 
equations are proposed. It is seen that the six modes 
are closely related to the low-lying excitations of the sys- 
tem. At finite temperatures, the magnetization, suscep- 
tibility and specific heat are calculated, which show var- 
ious behaviors for different couplings. The double-peak 
structure of the susceptibility and specific heat can be 
procured, whose positions and heights are found to be 
dependent on competing couplings. It has been shown 
that the XXZ anisotropy of F and AF couplings can have 
remarkable effect on the physical behaviors of the sys- 
tem. In addition, the experimental susceptibility, specific 
heat and magnetization of the diamond chain compound 
Cu3(C03)2(OH)2[6] can be nicely fitted by our TMRG 
results. 

For the spin- 1/2 frustrated Heisenberg diamond chains 
with AF couplings Ji, J2 and J3, the magnetization 
plateau at m = 1/6 in the ground state coincides with 
a perfect fixed sequence of the averaged local magnetic 
moment such as {..., {Sa, Sa, Sb), ■■■} with 2Sa + Sb — 1/2, 
which might be described by trimerized states. On the 
other hand, the static structure factor S{q) shows peaks 
at wave vectors q — 0, 7r/3 (57r/3), and 27r/3 (47r/3) 
for different external fields and different AF couplings. 
We note that the similar behavior of S{q) has been 
experimentally observed in diamond-typed compound 
Sr3Cu3(P04)4 In addition, the DMRG results of 

the zero-field static structure factor can be nicely fitted 
by a linear superposition of six modes. It is observed 
that the six modes are closely related to the low-lying 
excitations of the present case. At finite temperatures, 
the magnetization m{H), susceptibility x(^) ^^d specific 
heat C(T) demonstrate different behaviors at different 
AF couplings, say, the magnetization plateau at to = 1/6 
is observed whose width is found to be dependent on the 
couplings; the double peak structure is observed for the 
susceptibility x(T) and specific heat C(T) as a function of 
temperature, and the heights and positions of the peaks 
are found to be dependent on the AF couplings. 

For the spin-1/2 frustrated Heisenberg diamond chains 
with F couplings Ji, J3 and AF coupling J2, the mag- 
netization plateau at to = 1/6 in the ground state cor- 
responds to a perfect fixed sequence of the averaged lo- 
cal magnetic moment such as {..., {Sa, Sb, St), ■■■} with 
Sa + 2Sb = 1/2, which could be understood by trimerized 
states. The static structure factor S{q) shows peaks also 
at wave vectors q = 0, 7r/3 (57r/3), and 27r/3 (47r/3) for 
different external fields and different F couplings Ji, J3 
and AF coupling J2, which is expected to be experimen- 
tally observed in the related diamond-type compound. 
In addition, the DMRG results of the zero-field static 
structure factor can be nicely fitted by a linear superpo- 
sition of six modes with the fitting equations mentioned 
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above. The six modes are closely related to the low- 
lying excitations of the system. At finite temperatures, 
the magnetization m{H), susceptibility x{T) and spe- 
cific heat C{T) demonstrate various behaviors for differ- 
ent couplings, namely, the magnetization plateau at to = 
1/6 is observed whose width is found to depend on the 
couplings; the double-peak structure is also observed for 
the susceptibility x{T) and specific heat C{T), and the 
heights and positions of the peaks are found dependent 
on F couplings Ji, J3 and AF coupling J2. 

For the spin-1/2 non-frustrated Heisenberg diamond 
chains with AF couplings Ji, J2 and F coupling J3, the 
magnetization plateau at to = 1/6 in the ground state 
coincides with a perfect fixed sequence of the averaged 
local magnetic moment such as {..., (S'a, Si,, S'f,), ...} with 
Sa + '2Sb — 1/2, which could be understood by trimerized 
states. The static structure factor S{q) is observed to ex- 
hibit the peaks at wave vectors q = and 2tt/3 (47r/3) 
for different external fields and different AF couplings 
Ji, J2 and F coupling J3, which could be experimen- 
tally detected in the related diamond-type compound. In 
addition, it is found that the zero-field spin correlation 
function {SjS§) is similar to that of the = 1/2 Heisen- 
berg AF chain. At finite temperatures, the magnetiza- 
tion m{H), susceptibility x{T) and specific heat C{T) 
are found to reveal different behaviors for different cou- 
plings, i.e., the magnetization plateau at to = 1/6 is ob- 
tained, whose width is found to depend on the couplings; 
the double-peak structure is observed for the tempera- 
ture dependence of the susceptibility x(T) and specific 
heat C (T) , where the heights and positions of the peaks 
depend on different AF couplings Ji, J2 and F coupling 

The effect of the anisotropy of the AF and F inter- 
actions on the physical properties of the non-frustrated 
Heisenberg diamond chain is also investigated. For the 
case of the couplings satisfying Ji : J22 : ^3 = 1 : 
2 : —0.5, when the anisotropic ratio 72 = Jixj J2z = 
Jiyl J2Z 7^ 1, it is found that the width of the plateau at 
TO — 1/6, the saturation field, and the susceptibility x(T) 
show the same tendency, but quantitatively different, un- 
der the external field H along the z and x directions, 
while the specific heat C(r) for H along the z direction 
coincides with that along the x direction. For the case 
of the couplings satisfying J\ : Ji -Hz = 1 : 2 : — 0.5, 
when the anisotropic ratio 73 = Jzxl Jzz = Jzyl Jzz 7^ 1, 
it is seen that the width of the plateau at m = 1/6, the 
saturation field, and the susceptibility x(T) exhibit the 
opposite trends for H along the z and x directions, while 
the specific heat C(r) for H along the z direction also 
coincides with that along the x direction. 

For all the three cases, plateau states of m = 1/6 are 
observed during the magnetization, whose static struc- 
ture factor 5(5) shows peaks at wavevectors g = 0, 27r/3 
and 47r/3. But in absence of the magnetic field, the static 
structure factor S'(g) in the ground state displays peaks 
at g = 0, 7r/3, 27r/3, tt, 47r/3, and 57r/3 for the frustrated 
case with Ji, J2, J3 > 0; peaks at g = 0, 7r/3, tt, and 



57r/3 for the frustrated case with Ji, J3 < 0, J2 > 0; 
and a peak at g = tt for the non- frustrated case with Ji, 
'h > 0, J3 < 0. In addition, the DMRG results of the 
zero-field static structure factor can be nicely fitted by a 
linear superposition of six modes, where the fitting equa- 
tion is proposed. At finite temperatures, the double-peak 
structure of the susceptibility and specific heat against 
temperature can be obtained for all the three cases. It 
is found that the susceptibility shows ferrimagnetic char- 
acteristics for the two frustrated cases with some cou- 
plings, while no ferrimagnetic behaviors are observed for 
the non-frustrated case. 

The compound CU3 (003)2 (0H)2 is regarded as a 
model substance for the spin-1/2 Heisenberg diamond 
chain. The 1/3 magnetization plateau and the two broad 
peaks both in the magnetic susceptibility and the spe- 
cific heat have been observed experimentally!^. Our 
TMRG calculations with Ji : J2 : Jsz = I : 1.9 : -0.3 
and Jsx/Jsz = Jay/Jsz — 1-7 capture well the main 
characteristics of the experimental susceptibility, spe- 
cific heat and magnetization, indicating that the com- 
pound Cu3(C03)2(OH)2 may not be a spin frustrated 
magnet [26j. 



APPENDIX A: LOW-LYING EXCITATIONS OF 
SPIN-1/2 FRUSTRATED HEISENBERG 
DIAMOND CHAINS 

In this Appendix, the low-lying excitations of the spin- 
1/2 frustrated Heisenberg diamond chain are investigated 
by means of the Jordan- Wigner (JW) transformation. 
The Hamiltonian of the system reads 

N 

= ^^(•/lS3i-2 ■ ^3i-l + J2^3i-1 ' + J3^3i-2 ' ^3i 
i=l 

3N 

+ J3^3i-1 ■ ^3i+l + JlS3i ■ Ssi+i) — Sj,{Al) 

3 = i 

where 3iV is the total number of spins in the diamond 
chain, Jj > (i = 1,2,3) represent the AF couphng 
while Ji < the F interaction, and H is the external 
magnetic field along the z direction. In accordance with 
the spin configuration of the diamond chain, we start 
from the Jordan- Wigner ( JW) transformation with spin- 
less fermions 

j-i 

m— 1 

^1 = 4% - \. (A2) 

where j = 1, • • • , . Because the period of the present 
system is 3, three kinds of fermions in moment space can 
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be introduced through the Fourier transformations 

= ^E^''^''^«3fc. (A3) 

Ignoring the interactions between fermions, the Hamilto- 
nian takes the form of 

k 

+ (7iaifeajj, + 72a2fca|J^, + 7303^05*^, + /i.c.)](A4) 

where = f (2Ji + J2 + 2J3-6iI), wi=-(Ji + Ja)-^, 
^2=-i(Ji + J2+J3)-i?, ^3=^2, 7i = (Jie'"+J3e"'''')/2, 
72=(J2e'^)/2, and 73=(J3e*'^- + Jie-*'^')/2- 
Via the Bogohubov transformation 

aik uii(fc)aifc + ui2(fc)Q:2fc + ui3(fc)a3fc, 
a2fc = U2i(fc)aifc + -U22(fc)a2fc + U23(fc)a3fc, 
flsfc = U3i(fc)aife + W32(fc)a2fc + U33(fc)a3fc, (A5) 

the Hamiltonian can be diagonalized as 

3 

^ = ^9 + E E ^-^^tk^^^ ■ (A6) 

k i=\ 



The coefhcients of the Bogohubov transformation can be 
found through equations of motion ihhik = \aik^li\: 

^1 71 73 \ / Mli \ / Wii \ 

7l ^^2 72 Uli = Eifc W2i . (AT) 

73 72* W3 / \ M3j / \ W3, y 

For a given fc, the eigenvalues £ik and eigenvectors 
(wii, U2j, W3i) can be numerically calculated by the driver 
ZGEEV.f of the LAPACK, which is available on 
the website (29} . Figs [7] show the zero-field low-lying 
fermionic excitation e{k) for the frustrated diamond 
chain with different AF coupling, while Figs llSI present 
the zero-field low-lying fermionic excitation e{k) for the 
frustrated diamond chain with Ji, J3 < and J2 > 0. 
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